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We study the single field slow-roll infiation models that better agree with the available CMB and 

LSS data including the three years WMAP data: new inflation and hybrid inflation. We study 

these models as effective field theories in the Ginsburg-Landau context: a trinomial potential turns 

out to be a simple and well motivated model. The spectral index Us of the adiabatic fiuctuations, 

\^ , the ratio r of tensor to scalar fluctuations and the running index dus/dhik are studied in detail. 

(^ ■ We derive explicit formulae for Ua, r and dua/dlnk and provide relevant plots. In new inflation, 

c7^ ' and for the three years WMAP and 2dF central value n^ — 0.95, we predict 0.03 < r < 0.04 and 

CN I -0.00070 < dus/dlnk < -0.00055. In hybrid infiation, and for n^ = 0.95, we predict r ~ 0.2 and 

^r>. dus/dln k ~ —0.001 . Interestingly enough, we find that in new infiation ns is bounded from above 

r-{ ' by ris max ~ 0.961528 . . . and that r is a two valued function of Us in the interval 0.96 < Us < n, max. 

In the first branch we find r < rmax = 0.114769 .... In hybrid infiation we find a critical value fiQ ^^it 

for the mass parameter /io of the field a coupled to the infiaton. For fig < Ao Mp;/192, where Aq 

l/^ ' is the cosmological constant, hybrid infiation is ruled out by the WMAP three years data since it 

^S) ' yields a blue tilted Us > 1 behaviour. Hybrid inflation for fj,Q > Ao Mp;/192 fulIfiUs all the present 

, CMB+LSS data for a large enough initial infiaton amplitude. Even if chaotic infiation predicts Us 

CO ' values compatible with the data, chaotic infiation is disfavoured since it predicts a too high value 

^ ' r ~ 0.27 for the ratio of tensor to scalar fluctuations. The model which best agrees with the current 

^O , data and which best prepares the way to the expected data r < 0.1, is the trinomial potential with 

CO . negative mass term: new infiation. 
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I. INTRODUCTION AND RESULTS 

Inflation was introduced to solve several outstanding problems of the standard Big Bang model [l| and became an 
important part of the standard cosmology. At the same time, it provides a natural mechanism for the generation of 
scalar density fluctuations that seed large scale structure, thus explaining the origin of the temperature anisotropies in 
the cosmic microwave background (CMB), as well as that of tensor perturbations (primordial gravitational waves) [2,yl- 

A distinct aspect of inflationary perturbations is that these are generated by quantum fluctuations of the scalar 
field(s) that drive inflation. After their wavelength becomes larger than the Hubble radius, these fluctuations are 
amplifled and grow, becoming classical and decoupling from causal microphysical processes. Upon re-entering the 
horizon, during the matter era, these classical perturbations seed the inhomogeneities which generate structure upon 
gravitational collapse|2,|3J. A great diversity of inflationary models predict fairly generic features: a gaussian, nearly 
scale invariant spectrum of (mostly) adiabatic scalar and tensor primordial fluctuations, making the inflationary 
paradigm fairly robust. The gaussian, adiabatic and nearly scale invariant spectrum of primordial fluctuations provide 
an excellent fit to the highly precise wealth of data provided by the Wilkinson Microwave Anisotropy Probe (WMAP) [J, 
Q Perhaps the most striking validation of inflation as a mechanism for generating superhorizon ('acausal') fluctuations 
is the anticorrelation peak in the temperature-polarization (TE) angular power spectrum at I ^ 150 corresponding 
to superhorizon scales y. The conflrmation of many of the robust predictions of inflation by current high precision 
observations places inflationary cosmology on solid grounds. 

Amongst the wide variety of inflationary scenarios, single field slow roll models provide an appealing, simple and 
fairly generic description of infiation. Its simplest implementation is based on a scalar field (the inflaton) whose 
homogeneous expectation value drives the dynamics of the scale factor, plus small quantum fluctuations. The inflaton 
potential, is fairly flat during inflation. This flatness not only leads to a slowly varying Hubble parameter, hence 
ensuring a sufficient number of e- folds, but also provides an explanation for the gaussianity of the fluctuations as well 
as for the (almost) scale invariance of their power spectrum. A flat potential precludes large non-linearities in the 
dynamics of the fluctuations of the scalar field. 

The current WMAP data seem to validate the simpler one-field slow roll scenario |J,13- Furthermore, because the 
potential is flat the scalar field is almost massless, and modes cross the horizon with an amplitude proportional to 
the Hubble parameter. This fact combined with a slowly varying Hubble parameter yields an almost scale invariant 
primordial power spectrum. The slow- roll approximation has been recently cast as a 1/Nf, folds expansion |9(, where 
Nefoids ^ 50 is the number of efolds before the end of inflation when modes of cosmological relevance today first 
crossed the Hubble radius. 

The observational progress permit to start to discriminate among different inflationary models, placingstringent 
constraints on them. The upper bound on the ratio r of tensor to scalar fluctuations obtained by WMAP |jja| rules 
out the massless (f/^ model and necessarily implies the presence of a mass term in the inflaton potential [3, Q ■ 

Besides its simplicity, the trinomial potential is a physically well motivated potential for inflation in the grounds of 
the Ginsburg-Landau approach to effective field theories (see for example ref.^|). This potential is rich enough to 
describe the physics of infiation and accurately reproduce the WMAP data |j,l3- 

The slow-roll expansion plus the WMAP data constraints the inflaton potential to have the form Q 

F(0) = iVe/oWs A/4 «;(X) , (1.1) 

where (j) is the inflaton field, x is a dimensionless, slowly varying field 

(1.2) 



yNefolds Mpi 

w{x) ~ 0(1) and M is the energy scale of inflation which is determined by the amplitude of the scalar adiabatic 
fluctuations [J to be 

M ~ 0.00319 Mpi = 0.77 x lO^^GcV . 

Following the spirit of the Ginsburg-Landau theory of phase transitions, the simplest choice is a quartic trinomial for 
the inflaton potential [3, [l3 ■ 

f \ ,'^ 2 , It- [y 3 , y 4 /I ON 



where the coefficients wq, h and y are dimensionless and of order one and the signs ± correspond to large and small 
field inflation, respectively (chaotic and new inflation, respectively). Inserting ea. (|1.3|) in ea. l|l.l|l yields, 

V{^) = 1-0 ± — 02 _^ ^ 03 + - 04 . (1.4) 

where the mass m? and the couplings g and A are given by the following see-saw- like relations, 

U I ^ I \ \ ^ I \ T/ AT 1\/T^ 



where N = N^foids- Notice that y ~ 0{1) ^ h guarantee that g ^ 0(10^®) and A ~ 0(10^^^) without any flne tuning 
as stressed in ref. ^. That is, the smallness of the couplings directly follow from the form of the inflaton potential 
eq. ()l.l(l and the amplitude of the scalar fluctuations that fixes M Q . 

The small coupling Hmit y — > of eas. (|1.3() - (|1.4l) corresponds to a quadratic potential while the strong couphng 
limit y —f oo yields the massless quartic potential. The extreme asymmetric limit \h\ —^ (X) yields a massive model 
without quadratic term. In such limit the product \h\ M^ must be kept fixed since it is determined by the amplitude 
of the scalar fluctuations. 

We study here new inflation with the trinomial potential eas. (|1.3|l - H1.4|l and hybrid inflation [see below], the two 
models fulfill the observational constraints. We compute in both scenarios ng, r and the running dns/dlnk as 
functions of the parameters of the models, derive explicit formulae for Ug, r and dus/dlnk and provide relevant plots. 
Moreover, we plot the ratio r and the running dris/dlnk as functions of the scalar index rig. Since the value of n^ is 
now known Q-H, these plots allow us to predict the values of r and dus/dlnk for the different inflationary models 
considered. These predictions and plots are solely produced from theory and not from any fitting of the data. 

The three years WMAP data indicate a red tilted spectrum (n^ < 1) with a small ratio r < 0.28 of tensor to scalar 
fluctuations [5J . The present data do not permit to find the precise values neither of the ratio r nor of the running 
index drig/dlnk, only upper bounds are obtained |j,|a]- We therefore think that the value of Us [ea. (|4.2|l ] obtained 
through a fit of the data assuming r = dus/dlnk = is more precise than the values of Us obtained through fits 
allowing both r and dus/dlnk to vary. Notice that Us = 0.95 was independently found from the 2dF data under 
similar assumptions Q . More precisely, from the three years WMAP data p] as well as ref. ^ we take 

n^ = 0.95 ±0.02. (1.6) 

We find that for n^ — 0.95 and any value of the asymmetry h [see figs. 0]and|3], new inflation with the trinomial 
potential eas. (|1.3|l - H1.4() predicts 

trinomial potential new inflation for n,, = 0.95 : 0.03 < r < 0.04 and - 0.00070 < dn^/d In fc < -0.00055 . 
We find for the lower value n., = 0.93 of the three years WMAP data band: 

trinomial potential new infiation for 71^ = 0.93 : 0.003 < r < 0.015 and - 0.0011 < dn^/d In fc < -0.00033 . 

Moreover, in new inflation with the trinomial potential, we find that Hs is bounded from above by 

new inflation : Ug < Ug maximum = 0.961528 .... 

For Ug — 0.961528 ... we have in this model r ~ 0.114769 . . . (see figs. QJandEJl. Interestingly enough, there exists 
two values (two branches) of r for one value of Ug in the interval 0.96 < Ug < 0.961528. . . [see fig. 0]. The value 
Tmax — 0.114769 ... is the maximun r in the first branch. The values 0.16 > r > 0.114769 . . . correspond to a second 
branch of r as a function of Ug in the interval 0.96 < Ug < 0.961528 .... In the first branch we have 

r^a, =0.114769... . 

The absolute maximun value Vabg max ~ 0.16 belongs to the second branch and corresponds to the quadratic monomial 
potential obtained from ea. (|1.3() at y — 0. 

These predicted values of the ratio r fuUfil the three years WMAP bound including SDSS galaxy survey Q 

r < 0.28 (95% CL) . (1.7) 



Moreover, one can see from fig. 14 in ref. ^ that r < 0.1 (68% CL) from WMAP+SDSS. 

Chaotic inflation with the trinomial potential ea. (|l.^^() - (|1.4(l yields larger values of r than new inflation for a given 
value of Us Ji!]. More precisely, for Ug — 0.95 we flnd r = 0.27 for the binomial potential lOj (the trinomial potential 
introduces very small changes). 

Therefore, although the WMAP value for Ug [ea. ljl.GIl ] is compatible both with chaotic and new inflation, the 
WMAP bounds on r clearly disfavour chaotic inflation. New inflation easily fulfils the three years WMAP bounds 
on r and prepares the way for the expected data on the ratio of tensor/scalar fluctuations r < 0.1. 

In the inflationary models of hybrid type, the inflaton is coupled to another scalar field ctq with mass term — ^q < 
through a potential of the type [I3 

T/ fA„\ "^^ a2 , .9o ^2 ,2 , Mo f_2 4A 

Vhyb[(P,CFo) = 



2^2^ " I6A0 V M5 

= ^ c^' + Ao + ligl c^'-t^l) -i + J^^4 , (1.8) 

where m^ > 0, Aq > plays the role of a cosmological constant and gg couples ctq with (p. 

The initial conditions are chosen such that ctq and &o are very small (but not identically zero) and therefore infiation 
is driven by the cosmological constant Aq plus the initial value of the inflaton 0(0). The inflaton field 0(t) decreases 
with time while the scale factor a{t) grows exponentially with time. The field ctq has an effective classical mass square 



m 



2 2/2 ,,2 



3o r - Mo • (1-9) 



Since the infiaton field (p decreases with time, rn^ becomes negative at some moment during inflation. At such 
moment, spinodal (tachyonic) unstabilities appear and the field ctq starts to grow exponentially. Inflation stops when 
both fields (f> and ctq are comparable with cj) and &o and close to their vaccum values. 

We find that the time when the effective mass of the field ctq ea. H1.9|l becomes negative depends on the values of 
/Ug and 30 0^(0). For low values of ^g the field ctq starts to grow close to the end of inflation. On the contrary, for 
higher values of fi^ the field ctq starts to grow well before the end of inflation. This is explained by the fact that the 
scale of time variation of (Tq goes as ^q^. ctq evolves slowly for small /io and fastly for large /io [see flgs. 1711 1| . 

Only at A = hybrid inflation becomes chaotic inflation with the monomial potential ^ . For any value of 
A > even very small, the features of hybrid inflation remain. 

We compute Us-, r and dns/dlnk for hybrid inflation as functions of the parameters in the potential eg. p. 8(1 and 
the initial value of the inflaton field [see figs. I12I23| . 

The results of our extended numerical investigation of hybrid inflation can be better expressed in terms of the 
dimensionless variables 

A = ^^0 y = 'A ^^^ .2 ^ Mo M^l Nefolds 

M^N.foids ' ^Nefoids A Mpi ^ 2Ao ■ 

We depict in figs. I12l23l the observables ris, '' and the running index dns/dlnk as functions of A and Ug. We present 
a complete picture for hybrid inflation covering two different, blue tilted and red tilted, regimes. We find that for all 
the observables, the shape of the curves depends crucially on the mass parameter fi^ of the a field and the (rescaled) 
initial amplitude x{0) of the inflaton field. 

We find a blue tilted spectrum {ug > 1) for /i^ < /t^^j^ ~ 0.13 while for /i^ > fi^^^^ we can have either Ug > 1 or 
Ug < 1 depending on the initial conditions: for x(0) > x{0)crit we have Ug > 1, and for x(0) < x{0)crit we have Ug < 1. 
The value of x{0)crit grows with /i^: for fi'^ = 0.5, we find x{0)crit = 2.7 and for /i^ = 1.7, we find x{0)crit — 5.8. 

We see that Ug > 1 happens when the cosmological constant Aq is large enough compared with /Xg Mp^ N^foids- 
More precisely, p? < /t^^jj for Ag > 192 ^g Mpi using Nefoids = 50. That is, for Ag < 192 /Zg Mpi we have either red 
or blue tilted spectrum as explained above. 

For large A, Ug — \, r and dug/dlnk always tend asymptotically to zero whatever be ft^ and x(0). 

We see from our calculations that all blue tilted values of {ug, r) in the domain 1 < 71^ < 1.15, < r < 0.2 can be 
realized by the hybrid inflation model ea. (|1.8l) . However, at the light of the three years WMAP data ref. p| the blue 
tilted regime in hybrid inflation /z^ < ^J-'^j.^^ is ruled out. 

The situation is totally different in the red tilted regime /t^ > fL^^j^^ — 0.13 in hybrid inflation. The possible values 
of {ug , r) for such regime of hybrid inflation are in the upper-right quadrant as shown in fig. 1241 



Hybrid inflation in the red tilted regime fi^ > fJ.'^^u ^^'^ x(0) < x{^)crit fulfills the three years WMAP value for Ug 
[see ea. (|4.2|l ] as well as the bound on the ratio r [ea. l4.3|l ]. We can read from fig. 1201 and 1^ that 

/i^ > Aicrrf hybrid inflation : 0.2 > r > 0.14 and - O.OOK dn^/dln/s < for 0.952 < n^ < 0.97 . 

Notice that hybrid inflation in the red tilted regime yields a too large ratio r > 0.2 for Ug < 0.95. 

At the central three years WMAP value Ug = 0.95 both new and hybrid inflation are allowed. However, for Us < 0.95 
hybrid inflation is in trouble (r > 0.2) while for Us > 0.962 new inflation is excluded. 

The potential which best agree with the present red tilted spectrum and which best prepares the way to the expected 
data (a small r < 0.1) is the trinomial potential eas. (|1.3(l - (|1.4(l with negative mass term, that is small field (new) 
infiation. Hybrid infiation with a trinomial potential can also reproduce the present data in the red tilted regime 
fi^ > nl„t and x(0) < x{0)c„f 

All calculations presented in this paper stem from the inflaton potential in the slow roll approximation (dominant 
order in 1/A^ ~ 1/50). They do not use observational data as input. The analytical formulas and plots provided in the 
paper allow to read directly the predicted values of r and dris/dliik as functions of n^. In order to make illustrative 
predictions, we take the value Ug = 0.95 ± 0.02, as a judicious choice. The reader can see directly from the plots 
presented here our predictions for r and dus/dlnk for future observational values of n^. 

II. THE INFLATON POTENTIAL AND THE 1/Nefoids SLOW ROLL EXPANSION 

The description of cosmological inflation is based on an isotropic and homogeneous geometry, which assuming flat 
spatial sections is determined by the invariant distance 

(2.1) 



ds^ = dt^ 


- a^{t) dx^ . 


The scale factor obeys the Friedman equation 




1 da 

_a{t) dt 


' P{t) 
3 Ml, ' 



(2.2) 
a^tj at J o -"'-'pi 

where Mpi = l/%/8^ = 2.4 x lO^^ GeV. 

In single field inflation the energy density is dominated by a homogeneous scalar condensate, the inflaton, whose 
dynamics is described by an effective Lagrangian 



C = a^{t) 






(2.3) 



The infiaton potential V{(j)) is a slowly varying function of 4> in order to permit a slow-roll solution for the inflaton 
field <f,{t). 

We showed in ref. |3] that combining the WMAP data with the slow roll expansion yields an inflaton potential of 
the form 

T/(</.) = N M^ wix) , (2.4) 

where x is a dimensionless, slowly varying field 

X = -r=^ , (2.5) 

Vn Mpi 

w(x) ^0(1), 7V^50is the number of efolds since the cosmologically relevant modes exited the horizon till the end 
of inflation and M is the energy scale of infiation 

The dynamics of the rescaled field x exhibits the slow time evolution in terms of the stretched dimensionless time 
variable, 

f Af^ 

(2.6) 



Mpi VN 



The rescaled variables x ^^^ t change slowly with time. A large change in the field amplitude cj) results in a small 
change in the x amplitude, a change in ~ Mpi results in a x change ~ l/vN. The form of the potential, ea. H2.4|l 
and the rescaled dimensionless inflaton field eg. 1)2.511 and time variable r make manifest the slow-roll expansion as a 
consistent systematic expansion in powers of 1/A^ |9(]. 

The inflaton mass around the minimum is given by a see-saw formula 






2.45 X 10^^ GeV 



PI 



The Hubble parameter when the cosmologically relevant modes exit the horizon is given by 

H = Vn mH-^ 1.0 X 10^"* GeV = 4.1 TO, 

where we used that Ti. ^ 1. As a result, m <C M and H ^ Mpi. A Ginsburg-Landau realization of the inflationary 
potential that fits the amplitude of the CMB anisotropy remarkably well, reveals that the Hubble parameter, the 
inflaton mass and non-linear couplings are see-saw-like, namely powers of the ratio M/Mpi multiplied by further 
powers of 1/A^. Therefore, the smallness of the couplings is not a result of fine tuning but a natural consequence 
of the form of the potential and the validity of the effective field theory description and slow roll. The quantum 
expansion in loops is therefore a double expansion on (H/Mpi) and 1/A^. Notice that graviton corrections are also 
at least of order (H/Mpi) because the amplitude of tensor modes is of order H/Mpi . We showed that the form of 
the potential which fits the WMAP data and is consistent with slow roll eas. l|2. 411 - 12. 5|l implies the small values for 
the inflaton self-couplings Q. 

The equations of motion in terms of the dimensionless rescaled field x and the slow time variable r take the form. 



dx 






(2.7) 



The slow-roll approximation follows by neglecting the -h terms in eas. (|2.7|l . Both w{x) and h{T) are of order N^ for 
large N . Both equations make manifest the slow roll expansion as an expansion in 1/N . 

The number of e- folds N[x\ since the field x exits the horizon till the end of inflation (where x takes the value Xend) 
can be computed in close form from eqs. (|2.7(l in the slow-roll approximation (neglecting 1/A^ corrections) 



N[x\ 
N 



wjx) 
w'ix) 



dx <1 



(2.8) 



The amplitude of adiabatic scalar perturbations is expressed as 

M 



I kad\ j^2 7r2 



Kx) 



The spectral index Us , its running and the ratio of tensor to scalar fluctuations are expressed as 



(2.9) 



n, — 1 = 

N 
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N w{x) 
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(2.10) 



In eas. (|2.8|l - 12.10|l the field x is computed at horizon exiting. We choose A^[x] = A^ = 50. 
Since, w(x) and w'{x) are of order one, we find from ea. (|2.9ll 



M 



2^/3 7^ 

N 



A 



k ad I 



1.02 X 10" 



(2.11) 



where we used TV ~ 50 and the WMAP value for \A[^IJ = (4.67 ± 0.27) x lO'^ ^. This fixes the scale of infiation 
to be 

M ~ 3.19 X 10"^ MpL ^ 0.77 x 10^^ GeV . 

This value pinpoints the scale of the potential during inflation to be at the GUT scale suggesting a deep connection 
between inflation and the physics at the GUT scale in cosmological space-time. 

We see that |?is — 1| as well as the ratio r turn out to be of order l/N^foids- This nearly scale invariance is a natural 
property of inflation which is described by a quasi-de Sitter space-time geometry. This can be understood intuitively 
as follows: the geometry of the universe is scale invariant during de Sitter stage since the metric takes in conformal 
time the form 

ds' = J^ [idvf - {dxf] . 

Therefore, the primordial power generated is scale invariant except for the fact that inflation is not eternal and lasts 
for Nefoids- Hence, the primordial spectrum is scale invariant up to l/NefoUs corrections. The values n^ = 1, r = 
and dns/dlnk = correspond to a critical point as discussed in ref. 9]. This a gaussian fixed point around which the 
inflation model hovers in the renormalization group sense with an almost scale invariant spectrum of scalar fluctuations 
during the slow roll stage. 

The WMAP results favoured single inflaton models and among them new and hybrid inflation emerge to be prefer- 
able than chaotic inflation |10l |. 

We analyze in the subsequent sections new inflation and hybrid inflation in its simple physical realizations within 
the Ginzburg-Landau approach (the trinomial potential') \UJ^ . 

III. SPECTRAL INDEX Us, RATIO r AND RUNNING INDEX ^ FOR NEW INFLATION WITH THE 

TRINOMIAL potential" 

We consider here the trinomial potential investigated in ref.pXJl 

n0) = ^o-^0^ + ^0^ + ^0^ (3.1) 

where m^ > and g and A are dimensionless couplings. 

The corresponding dimensionless potential w{x) takes the form 

^ .2 h [y 3 , y 4,2 



where the quartic coupling y is dimensionless as well as the asymmetry parameter h. The couplings in ea. (|3.1(l and 
ea. (|3.2fl are related by, 



'-^^2N 


\MpiJ 
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8N 


( ^^ V 

\MpiJ 


the constant F{h) is related to Vq by 


Im- 


Vo 







(3.3) 



The constant F{h) ensures that w{x+) = w'{x+) = at the absolute minimum x — X+ — \/ (^ + l^'^l) "^^ the 
potential w(x)- Thus, inflation does not run eternally. F{h) is given by 

8 



F{h) = lh^ + 4h^ + l + ^\h\A^ , A=Vh^+l. 
o o 

The parameter h reflects how asymmetric is the potential. Notice that w(x) is invariant under the changes x ~^ 
—X, h ^r —h. Hence, we can restrict ourselves to a given sign for h. Without loss of generality, we choose h < and 
shall work with positive fields x- 
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FIG. 1: New Inflation. Us as a function of logy for the asymmetry of the potential \h\ = 0,1,3,7 and 16, y being the 
dimensionless quartic coupling. The y — > limiting value Us — 1— j^ = 0.96 is /i-independent and corresponds to the monomial 
potential i m . 
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FIG. 2: New Inflation. The running dus/dlnk as a function of logy for the asymmetry of the potential \h\ = 0, 1, 3, 7 and 16, 



y being the dimensionless quartic coupling. The y ^ limiting value 
monomial potential i m^ 0^. 



-0.0008 is fe-independent and corresponds to the 



Notice that y ^ 0(1) ^ h guarantee that g ~ 0(10 ^) and A ^ 0(10 ^^) without any fine tuning as stressed in 
ref. 0. 

New inflation is obtained by choosing the initial field x in the interval (0, x+). The infiaton x slowly rolls down the 
slope of the potential from its initial value till the absolute minimum of the potential x+ ■ 

Computing the number of efolds from eq. (|2.8I) , we find the field x s-t horizon crossing related to the parameters y 
and h. It is convenient to define the field variable z: 

_ y 2 

z=^X ■ 
We obtain by inserting ea. H3.2|l for w{x) into ea. (|2.8|) and setting N[x] — N, 



y = z-2h^-l-2\h\A + ^\h\ {\h\ + A - V^) 



Ifi 

--\h\{A + \h\)A'\og 



}_(^ , V^-\h\ 



2Fih)log[y^{A~\h\)] . (3.4) 



z turns to be a monotonically decreasing function of y: z decreases from z = z+ = (A + |/i|)^ till z = when y 
increases from y — till y — oo. When ^/z -^ y/zlf, y vanishes quadratically, 



2/^^^+ 2 (Vi-ViT)' + 0([Vl-ViT] 



We obtain in analogous way from eas. H2.9(l and (|2.10(l the spectral index, its running, the ratio r and the amplitude 
of adiabatic perturbations, 

^1_q]L z{z + 2h^z-lf y_ iz + Ah^z-l 

N [F{h)-2z+lhz^/^ + z^f N F{h)-2z+lhzm + z^ ' ^'^ 

dus 2^2 {z + 2hy/z-\){h+^^/I) 



'zy 



d^^k iV2 [^(/j) _ 2 z + I /l 23/2 + ^2] 

24 2 2 (z + 2/iV^-l)4 



2 



y z 



iV2 " \V(h\-')y^^ 



l4 



i^(/l)-2z+|;iz3/2 + z2] 



2 (3z + 4/iVz-l)(z + 2/iV^-l^^ 
^^ "^ ^ [F(/i) - 2 z + I /i z3/2 + z2 



16^ -(- + ^"^-^)' 2 , (3.7) 

^ [F(/l)-2z+|/lz3/2 + z2]2 



.(■S) |2_ JXL (^JLV [f(/t)-2z+|/tz3/2^^2 1 

12 7r2 VAfpJ 2/2 ^ (2_^2/iV^- 1)2 



\^tL? = TT7Z2 ( 77- ) ^ J... ,ol /. .^2 ^ • (3-8) 



IV. NEW INFLATION WITH THE TRINOMIAL POTENTIAL CONFRONTED TO THE THREE 

YEARS WMAP DATA 

We plot rig, its running and r in figs. m5]and|21as functions of logy for various values of the asymmetry of the 
potential /i, y being the dimensionless quartic coupling. Figs. ^and|31depict r and the running dns/dlnk as functions 
of Us for various values of the asymmetry h. 

We see that generically Ug < 1 and dug/dlnk < for new inflation for all values of the couplings. 

In new inflation we have the absolute upper bound 

Q 

new inflation : r < Tabs max = "^ = 0-16 , (4.1) 

which is attained by the quadratic monomial potential obtained from ea. (|3.2|) at y = 0. On the contrary, in chaotic 
inflation r is bounded as 

o ^a 

chaotic inflation : 0.16 = ^r < r < — - = 0.32 . 

N N 
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FIG. 3: New Inflation, r as a function of log j/ for the asymmetry of the potential |/i| = 0, 1, 3, 7 and 16, j/ being the dimensionless 
quartic coupling. The absolute maximun value r = -^ = 0.16 is reached for y — and all h and corresponds to the monomial 
potential eg. 15.211 . 
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FIG. 4: New Inflation, r as a function of Ua for the asymmetry of the potential \h\ — 0,0.15,0.4,0.7 and 20. For a given 
Us, r monotonically and slowly decreases with increasing \h\. r = rijis) is not too sensitive to h. The maximun value of 

rig is ,2™"^*'"'"" = 0.961528 . . . and the corresponding r is Vmax ~ 0.114769 The maximun value of r is Tabs max ~ 0.16 

and corresponds to the quadratic potential setting y = in eg. 13. 21 . For n^ = 0.95 (the three years WMAP value), we flnd 
0.03 < r < 0.04. 
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FIG. 5: New Inflation. The running dua/dlnk as a function of ns for the asymmetry of the potential \h\ = 0,0.15, 0.4, 0.7 and 20. 
The running turns out to be always negative in new inflation. For Us < 0.96, the running dus/dlnk decreases with increasing 
\h\. The opposite happens for Ua > 0.96. In the last case the dependence on h is weak. We find dns/dlnk = —0.00077 . . . 
at the branch point Us = 0.961 . . . for all values of \h\. The point n^ = 1 — -^ = 0.96, -jj^ = "^ ~ —0.0008 is reached 
for all values of h and corresponds to the monomial potential eg. 15. 21 . For n^ = 0.95 (the three years WMAP value), we find 
-0.00070 < dus/dlnk < -0.00055. 
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FIG. 6: New Inflation. Maxima of Us plotted vs. the asymmetry of the potential \h\. The limiting value for large \h\ is 
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This bound holds for all values the asymmetry parameter h. The lower and upper bounds for r are saturated by the 
quadratic and quartic monomials, respectively. 

We see from fig. ^that n^ exhibits a single maximun n^ maximum{h) as a function of the quartic coupling y for fixed 
asymmetry h. In fig. Elwe plot rig maximum{h) as a function of h. Us maximum{h) monotonically increases with \h\ and 
rapidly reaches its limiting value Us maximum = 0.961528 . . .. The corresponding value for r is r = 0.114769 . . .. Values 
ns > ns maximum — 0.961528 . . . cannot be described by new inflation with the trinomial potential eas. (|3.1( l- (|3.2|l . 

We see from fig. 0] and |31 that both r and the running dng/dhik are two- valued functions of n^ in the interval 
0.96 < Us < 0.961528 . . .. That is, for each Ug in this range there are two possible values for r and for the running 
dus/dlnk. Therefore, we can cover the whole range of values 0.96 < n^ < 0.961528 . . . choosing the lower branch for 
r. We find for this branch 

r<r„ax = 0.114769... . 

This maximun value Vmax is well below the absolute maximun in new inflation Tats max — 0.16 [ea. H4.1|l ] which belongs 
to the second branch. 

The plots of the ratio r and the running dns/dlnk as a function of Ug show that these quantities are not very 
sensitive to the asymmetry h for a given value of Us . 

The three years WMAP ^ data as well as ref. |a| yield for Ug the value (see also refs. 01 and Q) 

n^ = 0.95 ±0.02. (4.2) 

For Ug — 0.95 and any value of the asymmetry h [see flg. 0], new inflation with the trinomial potential eas. (|3.1|) - (|3.2|) 
yields 

new inflation: 0.03 < r < 0.04 and - 0.00070 < dn^/dhifc < -0.00055 . 

New inflation with the trinomial potential always yield n^ below the maximun value rig maximum = 0.961528 . . .. For 
Ug = 0.961528 ... we have in this model r = 0.114769 . . .. These values of the ratio r fuUfil the three years WMAP 
bound including SDSS galaxy survey [g] 

r < 0.28 (95% CL) . (4.3) 

Moreover, one can see from fig. 14 in ref. 5j that r < 0.1 (68% CL) from WMAP+SDSS. 

Chaotic inflation with the trinomial potential ea. (|3.1|) - (|3.2l) yields larger values of r than new inflation for a given 
value of Ug [l^l- More precisely, we find r = 0.27 for Ug = 0.95 for the binomial potential in chaotic inflation [l^ (the 
trinomial potential introduces very small changes). 

Therefore, although the WMAP value for Ug [ea. H4.2|l ] is compatible both with chaotic and new inflation, the 
WMAP bounds on r clearly disfavour chaotic inflation. New inflation easily fulfils the three years WMAP bounds 
on r. 

The present data do not permit to find the precise values neither of the ratio r nor of the running index dng/dln fc; 
only upper bounds are obtained j^i^]. We therefore think that the value of n^ [ea. (|4.2|) ] obtained through a fit of the 
data assuming r = is more precise than the values of Ug obtained through fits allowing both r and dug/dlnk to 
vary. Notice that Ug — 0.95 was independently found from the 2dF data under similar assumptions [3. 

Ref. I3] reports fits yielding negative values for dug/dln k of the order ^ —0.05. Notice that the order of magnitude 
of the running dug/dlnk is just fixed by the fact that it is a second order quantity in slow-roll: ^ ^ ^ 0.0004. Still, 
the negative sign of the running reported by ref.[5J agrees with the sign prediction of new inflation with the trinomial 
potential [see fig. [21and|S]. 

In summary, new infiation with the trinomial potential ea. H3.1(l - (j3.2(l predicts 0.03 < r < 0.04 and —0.00070 < 
dug/dlnk < -0.00055 for Ug = 0.95. For 0.93 < Ug < 0.962 it predicts 0.01 <r < 0.115 and -0.001 < dug/dlnk < 
-0.0003 [see figs. Eland[5]. 

V. LIMITING CASES OF THE TRINOMIAL POTENTIAL IN NEW INFLATION 

Let us now consider the limiting cases: the shallow limit {y — > 0), the steep limit y — > 00 and the extremely 
asymmetric limit \h\ — > 00 of the trinomial potential for new inflation eas. H3.1| l- H3.2(l . 
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A. The shallow limit y ^ of the Trinomial Potential 

In the shallow limit y —> 0, z tends to z = z+ = (A + \h\)'^, which is the minimum of y in ea. H3.4|l . We find from 
eqs.(|23l-(E3J), 

which coincide with rig, ^,"°^ and r for the monomial quadratic potential. That is, the y — + limit is /i-independent 
except for | A^ j^|. For fixed h and y ^ the inflaton potential ea. (|3.2|) becomes purely quadratic: 

w{x) '=° A(A + \h\) ix ~ x+f + 0{^) , (5.2) 

where x+ = \ - (^ + I^D- Notice that the amplitude of scalar adiabatic fluctuations ea. (|5.1|l turns out to be 

proportional to the square mass of the inflaton in this regime which we read from ea. (|5.2() : 2 A(A + \h\). The shift 
of the inflaton field by x+ has no observable consequences. 

The numerical values in ea. H5.1|l are in agreement with figs. 1 1131 in the y — > limit. For h = we recover the results 
of the monomial potential. 

B. The steep limit y ^ oo of the Trinomial Potential 

In the steep limit y —> oo, z tends to zero for new inflation. We find from ea. (|3.4fl 



y '^" -F{h) logz - q{h) - 1 + 0(yi) , (5.3) 



where 



q{h) = 2 F{h) log[A - \h\] -^ {h^ + \h\ A) <{ 8 A2 log 



2 I A 



q{h) is a monotonically increasing function of the asymmetry \h\ : < q{h) < oo for < \h\ < oo. 
Then, eqs.(E2I)-(|n7I) yield, 

y^i _ y y>i IGy _.«±l+2lM 

""' ^ NFih) ' "■ " NF^h)^ 

dris «>i 2 y"^ \h\ _y+i+£W 



dlnfc N^ F^{h) 



(5.4) 



In the /i — > limit we recover from eas. (|5.3l l- (|5.4ll the results for new inflation with a purely quartic potential: we 
have F(0) — 1 and g(0) = and ea. H5.4|l becomes, 

jj»i,/i^o _, y a>iJi^o 16 y ^ ^^ 
n, = 1 , r — e ^ , 

iV ' TV ' 

dus y»i,/i^o 2 2/2 \h\ _ 1 
— e " , 



dink N^ 

(S) |2V»M.^o N^ f M y e^ 
12^2 {mpJ y 






The behaviour in eas. (|5.4|) is in agreement with figs. Ill3l in the y — > +00 limit. 
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C. The extremely asymmetric limit \h\ ^ cx) of the Trinomial Potential 

Eas. H3.4| )- H3.8() have a finite limit for \h\ -^ oo with y and z scahng as h^. Define, 

z- ^ Y- y 

Then, we find for \h\ — > oo from eas. (|3.4|) - (|3.8|l keeping Z and Y fixed, 

4 r- 4 Z 16 

Y = Z -- VZ -4-- log- 



3 3^4 3\/Z ' 

y Z^{^-2f Y 3Z-4x/Z 

TV [M „ I ^f + z2]2 + iv f _ I zi + Z2 ' ^ ■ ^ 



dn, 2 ^, ^ iVz - 2){l VZ - 1) 



dink N^ r 16 „ 8 2-3/2 ^ ^21 2 

^V^Z^ (^-^)^' 8 ^,^.(3VZ-4)(VZ^2)2 



-/V^ [16 _ 8 r^3/2 I 721'' iV^ ri6 



[^ _ 8 ZV2 + ^2] * 7V2 [^ _ 8 ^3/2 + ^2] 



^Y Z2(VZ-2)2 
r = 16 — 



iv [f^izi + z^Y 



' '=-'^' 12^2 yMpJ Y^Z^{VZ-2)^ ■ ^ ^ ^ 

We have < Z < 4 for +oo > F > 0. In the |ft,| — > oo fimit the inflaton potential takes the form 

w{x) _ 32 1 ^/r ^^3 y 4 



^(^) = |J™oo^^ = 3y-3VT^ +32^ 

This is a broken symmetric potential without quadratic term. Notice that the cubic coupling has dimension of a mass 
in ea. H3.1|l and hence this is not a massless potential contrary to the quartic monomial x^- In addition, ea. H5.8|) shows 
that for large \h\ one must keep the product \h\ Af2 fixed since it is determined by the amplitude of the adiabatic 
perturbations. We see from ea. H5.8|l that AI = \/\h\ M becomes the energy scale of inflation in the \h\ -^ oo limit. 
M ~ lO^^GeV according to the observed value of |A^ ad\/^ displayed in ea. (|2.11|) . while M and m, 

M.^"'l=-0 , .. = ill = ^ 1"'^- . 

yn Mpi \h\ Mpi 



vanish as \h\ — j oo. 

The curves in figs. Ill6l for high values of |/i| are well described by ea. H5.6|l . 



VI. HYBRID INFLATION 



In the inflationary models of hybrid type, the inflaton is coupled to another scalar field ctq with mass term — /Iq < 



through a potential of the type 



lyb r: 

m 



^^ ,2 , 90 J.2 2 , A^O f 2 4Ao^^ 



= ^ c^' + Ao + lial c^'-f^l) -'o + j^^'^'o , (6.1) 

where m^ > 0, Aq > plays the role of a cosmological constant and gg couples co with 0. 

The initial conditions are chosen such that co and ao are very small (but not identically zero) and therefore one 
can consider initially, 

2 

Vhyb{cj),0) = '^^^ + Ao. (6.2) 
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One has then inflation driven by the cosmological constant Aq plus the initial value of the inflaton 0(0). The inflaton 
field (/)(i) decreases with time while the scale factor a{t) grows exponentially with time. We see from ea. H().l|l that 



2 i2 

9o </> 



/^o 



(6.3) 



plays the role of a effective classical mass square for the field (Tq. The initial value of m^ depends on the initial 
conditions but is typically positive. In any case, since the inflaton field if) decreases with time, m^ will be necessarily 
negative at some moment during inflation. At such moment, spinodal (tachyonic) unstabilities appear and the fleld 
ctq starts to grow exponentially. Inflation stops when both fields (f> and ctq reach their vaccum values. A matter 
dominated regime follows. 

Normally, the field ctq is negligible when the relevant cosmological scales cross out the horizon. Hence, ctq does not 
affect the spectrum of density and tensor fluctuations except through the number of efolds. Hence, hybrid inflation 
is a single-field inflationary model as long as </> solely contributes to the spectrum of density and tensor fluctuations. 
However, we find specific regions of parameters {go, fj,o, Aq) and initial conditions where both fields (j) and ctq contribute 
to the cosmologically relevant fluctuations making hybrid inflation a two-field inflationary model. We do not consider 
such regions of parameters here which are outside the scope of this paper. 

In terms of the dimensionless fields and couplings, the potential Vhyb{4>, ctq) ea. (|().l|l reads 

,4 / n A X 2 



, , 1 2^ mV 2 2 A 



1 2 2 2 

25 ^ X 



ix^ + i(5^X^-.V^ + iA+|^.^ 



(6.4) 



where 



a{T) 



<Ja{t) 



o_.N Mil 



N M, 



3o 



,2^i^ ^^^ ^_ 



PI 



2Ao 

M'^ N 



The evolution equations for this potential in dimensionless variables take the form 



'^^H 


'-.^'-^'-f. 


( 2 2A 

"' 2 


1,222 
1 +.9 cr X 


\ d d TO 


x = o, 


N dr^^^ 


n±-,'..^. 


2A 


a(r) - . 



(6.5) 



Since the field a is chosen initially very small, it can be neglected and we can approximate the evolution equations 



(|^ as 



3WX + X-0 , 'n\T)=^- [x'+A] 



The number of efolds from the time r till the end of infiation is then given by ea. (|2.8|l . 



N{t) 



N 



n{T) dr 



xe„. ^(^) ^ r 2. ^ 2 1 _^ ^ , , X{r) 
—n—^ dx^ -r[x (r) - Xendl + ^ A log 

X(r) W'[X) 4 ■- -'2 Xend 



(6.6) 



(6.7) 



Xend is the inflaton field at the end of inflation. 

We see that the inflaton fleld and its dynamics only appear in N{t) ea. H6.7|l through the value of Xend where 
infiation stops. The value of Xend follows by solving eas. (|6.5|) and depends on the initial conditions as well as on the 
parameters g, fi and A. 

The spectral indices are given by eas. H2.1U|l and the amplitude of adiabatic perturbations by ea. l|2.9|) . By using the 
potential ea. H6.2|l in dimensionless variables we find. 



w'(x,0) = -(x'+A) 
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l^fcadl 9Qn^\MpiJ x^ ' iV(x2 + A)2' ^"■''^ 

_ 4 A-2x^ rfn, _ 32 x^(2A-x^) , . 

iV (x2 + A)2 ' dink N^ (x^+A)4 ' "^ ' ' 

where x is the inflaton at the moment of the first horizon crossing. 

VII. INFLATON DYNAMICS IN HYBRID INFLATION 

In figs. 1711 II we display the numerical solution of tlie equations of motion H().5|l as functions of time for A = 4 A^ = 
200 , /i^ — 1.7 A, 0(0) = 2.3 VA. We see first a stage of slow-roll quasi-de Sitter inflation tiU r ~ 39 in this example. 
Namely, \(j)\ <C |0| and h are practically constant during this lapse. In this slow-roll stage the equations of motion 
(|6.6() can be integrated in close form with the solution 



-^(r-ro) = Argtanh^i==- Jl + ^, (7.1) 

3A y^T^ V A 

which defines x = xi^) ^tnd where tq is an integration constant. Notice that x{^) is a monotonically decreasing 
function of time since x = — x/[3 h] < [ea. (|6.fi|l ]. 
When X ^ A, ea. H7.1|l approximates by 



x(t) ^Xoe V 3^ ^ , X < A , 
while in the opposite limit x 3> A from ea. (|7.1() we have, 

xM-Xi-v/gT, X>A • 

Here, xo f^nd xi are integration constants. 

We have verified that ea. (|7.1(l as well as ea. (|6.7(l provide an excellent approximation to the numerical solution of 
eqs. (|6.5|) . 

The number of efolds during inflation is about 280 in the example depicted in figs. I7I11I larger than the required 
minimun of about 60 efolds. This stage is followed by a matter dominated era. We choose a very small initial 
amplitude for the sigma field and its time derivative. The sigma field stays very small till its effective mass square 
[ea. l|6.3|l ] becomes negative and spinodal unstabilities show up. At this moment, (r ~ 4 in this example) the sigma 

4 

field as well as its time derivative start to increase exponentially fast till the growth of the non-linear term +yt "'^ 
in the last equation in ea. H6.5|l shuts off the unstabilities. 

Inflation stops at the moment when both x and a are comparable with x and ct (r ~ 39 in this example). At 
this time, both x and a are very close to their vaccum values Xvac — and <Ti,ac — ■ That is, when the kinetic 

terms become relevant, the energy is no more dominated by the vaccum energy. At the same time, the slow roll 
approximation ceases to be valid. 

The time when the effective mass of the field a [see eq. ()6.3|l ] 

ml^m^g'x'-^^') (7.2) 

becomes negative and a starts to grow depends on the values of /i^ and g^ X^(0)- For low values of fi^: (typically 
for /^^ < 0.08 A when x^(0) < A, and ^^ < 0.2 A when x^(0) < 2. A), the field a starts to grow close to the end of 
inflation. On the contrary, for higher values of /x^ (typically for /i^ > 0.08 A when x^(0) < A, and /i^ > 0.23 A when 
X^(0) < 2. A), the field a starts to grow well before the end of inflation. This is explained by the fact that the scale 
of time variation of a goes as /i~^; a evolves slowly for small /i and fastly for large /i. 

& exhibits a peak around the point where m^ changes sign and then returns to a very small value while a slowly 
approaches its vaccum value. This evolution is depicted in fig. 1101 

In the example depicted in figs. 1711 II the effective mass square [ea. H7.2|l ] of the a field changes sign at r '--^ 4 well 
before x reaches its vaccum value (zero) and infiation ends. This follows from the choice of a large value for ^^ in 
figs. 1711 II For smaller values of fi^, rrv^ [ea. (|7.2|) ] flips its sign later when the inflaton x is much smaller. 
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FIG. 7: Hybrid inflation. The logarithm of the scale factor (the number of efolds) as a function of time. The chosen parameters 
in ea. (l531l are A = 4 iV = 200 , n"^ = 1.7 A, 4>{0) = 2.3 \/A. A stage of slow-roll quasi de Sitter inflation takes place (tiU r ~ 39 
in this example) followed by a matter dominated era. 



70 



60 



50 



40 



30 



20 



10 





1 1 


1 


1 






1 




1 


1 




/ 
/ 


- 




















/ 


- 


- 
















/' 


/ 




- 




h vs. time — 












/ 


/ 
/ 






' 


- 


1/h vs. time --- 










/ 










- 


- 










r^' 












- 


- 


1 , 1- - 


^.^ 


1 

1 






1 




1 


1 




- 







10 



20 



30 



40 



50 



60 



70 



80 



FIG. 8: Hybrid inflation. The Hubble parameter h and its inverse 1/h as a function of time. Same parameters as in flg. 
h slowly decreases with time in the slow-roll quasi de Sitter stage (till r ~ 39 in this example) followed by /i ~ 2/ [3 r] in the 
matter dominated era. 
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FIG. 9; Hybrid inflation. The inflaton field x s-nd its time derivative as a function of time. Same parameters as in fig. |7| 
1x1 ^ 1x1 during tlie slow-roll inflationary stage. Inflation stops when X ~ X "^ (^t r ~ 39 in this example). 
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FIG. 10: Hybrid inflation. The field sigma a and its time derivative as a function of time for hybrid inflation. Same parameters 
as in fig. Q The fields a and a start with small values and grow exponentially fast when m^ < [ea. H7.2l ] (at r ~ 4 in this 
example). 
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FIG. 11: Hybrid inflation. The equation of state pressure/energy density as a function of time. Same parameters as in fig. 13 
The equation of state clearly shows the two stages: pressure/energy = — 1 during inflation followed by oscillations with zero 
average pressure in the matter dominated era. 
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FIG. 12: Hybrid inflation. The ratio r vs. n^ for fi' = 0.05 A < fi^^n , 5^ = j and 0.01 VA < x(0) < 0.4 VA. Notice that here 
r has a maximun as a function of ris- In addition, Ug > 1 for all values of A and x(0) since /i^ < ficrit 
with x(0) for this range of x(0)- 
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In order to compute the observables Ug, r and dus/dlnk from ea. H6.8|l we need the value of the inflation field x 
at 50 efolds before the end of inflation. We thus integrated numerically eas. (l6.5|) till the end of inflation and then 
extracted the value of x at 50 efolds before. We define the end of inflation as the point where the ratio pressure over 
energy reaches 10%. This gives t = Tend — 34 for the example in fig. ^] 

At A = hybrid inflation becomes chaotic inflation with the monomial potential ^ X^- We want to stress that only 
at A = hybrid inflation becomes chaotic inflation. For any value of A > (even very small) the features of hybrid 
inflation remain. The time Tend gets longer and longer for A ^ 0^. 

VIII. SPECTRAL INDEX n,, RATIO r AND RUNNING INDEX -^ IN HYBRID INFLATION 

* ' din k 

We see from eas. H6.9ll that the field x naturally scales as ^/A. It is then convenient to introduce the rescaled field 
and the rescaled mass 

X=^ ' ^ -A- (^-^^ 

Then, eas. l|(i.9|l take the form 

A 



^"(x,o) = ^(x' + i) 

<S) .,_N'A'fM V (X' + 1)=^ _ 32 x' 



^^kad\ - gQ^2 [mpJ x" ' ''^ NA ix^ + l)^' ^^-^^ 

r, -1 , 4 l-2x^ dn, ^ 32 x^(2-x^) 

7VA (x2 + l)2 ' dink 7V2A2 (x2 + l)4 • ^ ' 

Notice that (n^ — 1) may have either sign according to ea. (|8.3l) . Hybrid inflation is usually associated with red tilted 
spectrum (n<, > 1). However, both regimes, Us > I and tIs < 1 are realized by hybrid inflation. 

Whether n^ > 1 or n^ < 1 for a given set of parameters A, /x^, g and the initial conditions is a dynamical question 
that can only be answered after evolving the fields till the end of inflation according to eas. H6.5|l . As we see in ea. (|6.9|) 
the question is whether twice x^ at horizon exit is larger or smaller than A. Recall that horizon exit happens about 
50 efolds before the end of inflation and that x^ monotonically decreases during inflation. Even if initially 2 x^ > A, 
it can be very well that 2 x^ < A at horizon exit. This depends on how many total efolds Nt > 60 wc have; since 
horizon exit happens {Nt — 50) > 10 efolds after the beginning of inflation, the larger is {Nt — 50) the smaller can 
be x^ at horizon exit. 

We vary the parameters A, /z^ and the initial conditions always keeping the total number of efolds Nt during 
inflation larger or equal to 60. We keep g^ — j since this parameter is less relevant than the others. We explored the 
parameters region where r < 0.2 and 0.95 < n^ < 1.15. 

A. Red tilted and blue tilted regimes in Hybrid Inflation 

Extended numerical investigation showed that there exists a critical value of fj?, fJ-erit — 0-13 A such that Ug > 1 
provided fi"^ < ^^erit ■ 

For yu^ > Herit '^6 ^^'^ both regimes, rig > 1 and Ug < 1. This property is valid for all initial values of the inflaton 
compatible with the restrictions Nt > 60 and one-inflaton fluctuations. Otherwise, if the held a is relevant at horizon 
exit we should also include its contribution to the density fluctuations. Such calculation is beyond the scope of the 
present work where we concentrate on single inflaton fluctuations. 

The larger is ^^, the earlier inflation ends, the earlier horizon exit happens and the large is x at horizon exit. That 
is, increasing /i^ decreases Ug. This explains why we necessarily find rt^ < 1 for /i^ > ulrit- 

For fj,^ > ^erit "^^ fin.'i that Ug > 1 for x(0) > x{^)crit- That is, increasing x(0), increases Ug. This is so because the 
larger is x(0); the larger is Nt since Nt ~ x(0)^ [see eas. H6.6|l and H6.7|l ]. Then, the larger is {Nt — 50) the smaller 
is X at horizon exit and the larger is Ug. 

In all cases, (both fp > /i^^^j and fi^ < fJ-l^it) for A ^ oo we always find n^. ^ 1, r -^ and -j^ -^ 0. 
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FIG. 13: Hybrid inflation. The ratio r vs. Us for ^^ = 0.05 A < fi^^.t , 5^ = | and 0.5 \/A < x(0) < 2.4 VA. Notice that here 
r has a maximun as a function of n^. In addition, n^ > 1 for all values of A and x(0) since fi^ < ficrit — 0.13 A. r decreases 
with x(0) for this range of x(0). 
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FIG. 14: Hybrid inflation. The ratio r vs. Us for ^^ = 0.13 A ~ fi^^.t , 5^ = i and 0.7 VA < x(0) < 2.3 VA. Notice that here 
Us > 1 for all values of A and x(0). 



22 



0.014 



0.012 



0.01 



0.008 



0.006 



0.004 



0.002 



X(0) = 0.7 
X(0) = 0.9 
X(0) = 1.1 
X(0) = 1.3 
X(0) = 1.5 
X(0) = 1.7 
X(0) = 1.9 
X(0) = 2.1 
X(0) = 2.3 



1.02 



1.04 



1.06 



1.08 



1.1 



1.12 



1.14 



1.16 



FIG. 15: Hybrid inflation. The running dus/dlnk vs. Us for fi^ = 0.13 A ~ f4,.it , g^ = \ and 0.7 VA < x(0) < 2.3 \/A. 
Notice that here dUs/dlnk > for all values of A and x(0)- It exhibits a shape similar to r vs. n^ [see fig. 1131 . 
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FIG. 16: Hybrid inflation. The index Us vs. A for /i^ = 1.7 A > /i^.^, 5^ = 3 and 2.7 VA < x(0) < 5.7 \/A ~ xWcrit- Notii 
that here Us < 1 for all values of A and this range of x(0). 
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Figs. I12I23I show the observables n^, r and the running index dus/dhik as functions of A and n^ for /i^ = 
0.05 A, ^^ = 0.13 A and ^^ = 1.7 A. A complete picture for hybrid inflation emerges covering two different, blue 
tilted and red tilted, regimes. We find that for all the observables, the shape of the curves depends crucially on the 
mass parameter jl^ of the a field and the (rescaled) initial amplitude x{0) of the inflaton field. 

We find three regimes according to the value of jl^ : 



• 



fi^ < 0.075. Here we always have n^ > 1 and r has one maximun as a function of n^ (or A). Ug monotonically 
decreases with A. Figs. 1121 and [TSl show r vs. n^. for /i^ = 0.05 A and various values of x(0)- ^ displays a 
maximun as a function of n^. In addition, r grows with x(0) for x(0) < 0.5 while it decreases with x(0) for 
x(0) > 0.5. The running -j^^ behaviour is qualitatively similar to the behaviour of r above described. The 
running -j^fj: is here positive and grows when Ug grows. 



r and the running -j^^ vs. n^, respectively, for /l^ = 0.13 ~ Acrit- ^^ ^^^ ihaX the running -j^^ behaviour is 
qualitatively similar to the one of r. 



• 0.075 < jl < jj.'^rit — 0-13. Here we always have Ug > 1 and r monotonically grows with A. Ug monotonically 
decreases with A for fi^ < 0.1 while it exhibits a maximun as a function of A for p,^ > 0.1. Figs. ll4land[TKldcDict 

' and the runn 
qualitatively sii 

• /t^ > 0.13 ~ /t^^j^. Here rig > 1 for x(0) > x{0)crit and iig < 1 for x(0) < x(0)c«t- 

The value of x(0)crjt grows with fi^: for ji^ = 0.5, wc find x(0)crit = 2.7 and for /t^ = 1.7, we find x(0)crjt = 5.8. 

For x(0) < xiO)crit, ng monotonically increases with A with values Ug < 1. For x(0) > x(0)crit, ng shows an 
absolute mELximum, which is always Ug max > 1- The highest Ug values concentrate and narrow in the small 
A region. It must be noticed that for each curve, [each x(0)]j '^s can take values Ug > I and rig < 1: even if 
ng max > 1, 'T's can be below unit in the two sides of the curve, [see figs. 1161 and [T7 ] . 

The value Ug = 1, is reached asymptotically for large A from rig < I for both x(0) < x(0)crit and x(0) > x(0)crit- 
In addition, the value ng — 1 with 0.2 > r > 0.04 is found for a variety of values of A and x(0) > x(0)crit as we 
see from figs. [T71 [Bl and IHl 

B. The ratio r in Hybrid Inflation 

The ratio r in figs. 1181 and [T^ exhibits an oscillatory pattern and two different regimes: 

For x(0) < x(0)crit, (fig- ITH|l r decreases monotonically reaching very small values for large A. For x(0) < x{0)crit, 
r does not feature any oscillation. 

For x(0) > x{^)crit, r decreases with A, (fig. I19|) r has an absolute minimum r„iin and then grows till a maLximum, 
Tmax- The oscillations show up and concentrate with growing amplitude for small A for high x{Q), Tmin and r^ax shift 
towards the smaller A with increasing x(0); '"'min decreases, and r„iax increases, for increasing x(0)- The convexity of 
the curve for small A increases for decreasing x(0)- 

For x(0) > x(0)crit, each curve [each x(0)] shows for r a oscillatory behavior with three clear parts: (1) the 
asymptotic part of monotonically decreasing r for large A at the right of Tmax] (2) the increasing part at the left of 
Tm.ax ; (3) the sharp decreasing part for small A at the left of r„ji„ . In the minima, rmin can be extremely small for 
small A, which is a new feature in hybrid inflation. 

In the asymptotic regime of large A, r does not feature any oscillation. All curves [for all x(0)] coalesce into r = 
for A — > c». 

There are three distinct regimes: small A, intermediate A and large A. The new oscillatory behavior for high and 
intermediate x(0) > x(0)cT-it is in the region of small and intermediate A. The monotonically decreasing behavior for 
low x(0) < x{0)crit is in the asymptotic region of large A. 

The highest values of r appear for small A whatever be the hybrid regime; for such high values of r both regimes 
x(0) < x{^)crit and x(0) > x(0)crjt superpose. From such high values, r decreases sharply till its minimun rmm h^ 
the small A region for x(0) > x{^)criu or r decreases monotonically reaching asymptotically the large A regime for 
x(0) < x(0)crjt- The low values of r for low A are a totally new feature in hybrid inflation. 

C. The running index dns/d\nk in Hybrid Inflation 

The curves of the running index dug/dhik figs. 1151 1221 and 1231 show new features in two different regimes. 

For x(0) > x(0)crjt, the running index dug/dlnk shows a similar shape as r and is essentially positive. It oscillates 
with at least one maximum and one or two minima and three different components: (1) the asymptotic part of 
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monotonically decreasing running with increasing A, at the right of the maximum, going to zero in this regime; 
(2) the increasing running with A, which is a new feature in hybrid inflation, and (3) the sharp decreasing of 
the running tifl its minimum value for small A. The highest running appears for small A as in the known hybrid 
regime. The lower running values for small A, as well as the oscillations for small and intermediate A, are totally new. 

For low x(0) < x(0)crit the running index does not exhibit any oscillation. Both for x(0) < x{0)crit and x(0) > 
x(fi)crit dus/dlnk grows with A till it reaches its maximun and then decreases monotonically with A. dng/dlnk 
vanishes asymptotically for large A without any oscillation. This is a totally new feature for hybrid inflation. 

Thus, hybrid inflation describes both dus/dlnk > and dus/dlnk < 0. It must be noticed that dris/dhik < can 
reach very low values for small A which is a totally new feature in hybrid inflation. 

In summary, the new^ features for the running index in hybrid inflation are: both positive and negative running, 
increasing of the running with A, transition from positive to negative running passing through zero running when A 
grows. 

D. r vs. Us. Confrontation of Hybrid Inflation to the three years WMAP data 

r vs. Us in figs. 12012 II depicts a oscillatory behavior clearly showing the two regimes: x{0) > x{^)crit corresponding 
mainly to Ug > 1, although it also covers a small portion of ng < 1, and x(0) < x(0)crit for which Us is entirely red 
tilted. 

All curves end [x(0) < x{0)crit], or start [x(0) > x{0)crit], at rig — 1, which is the fixed point for all x(0)j with 
three different behaviors: 

(1) the sharp decreasing of r in the range Us < 1, approaching n^ = 1 as the end point, this is for x{0) < x{0)crit, 
in which r can take high values for small n^, (n^ near 0.95). 

(2) the monotonically decreasing of r with n^ at the right of r„iax, for x(0) > x(0)crit, in which r vanishes 
asymptotically for 'high' Ug , {ug > 1.07). 

(3) The new hybrid behavior for high x(0) > x{0)crit in between the above two regimes, in which r shows a 
maximum and a minimum between two sharp decreasing and increasing 'arms', lying at n^ < 1 and n^ > 1 respectively. 
rmin decreases and rmax increases as increasing x(0)- AH Vmax he in the red tilted regime n^ < 1. All rmin he in the 
blue tilted regime n^ > 1. 

For x(0) > x{0)crit, all curves go towards rig = 1, r = 0. Most of each curve lies in the Ug > 1 region, r 
monotonically decreases with Ug in the range Ug < 1 from the maximun r^nax going towards r = for ng —> 1. For 
x(0) < x{0)crit the curves pill up in the n^ < 1 region with almost the same slope, r sharply decreases in this region. 

The curves for x{0) < x{0)crit piU up in the n,, < 1 region with almost the same slope, r sharply decreases in this 
region. 

All the blue tilted values of {ng,r) in the domain 1 < n^ < 1.15, < r < 0.2 are realized by hybrid inflation. 

The red tilted regime in hybrid inflation can only be realized for fj? > li^^it — 0-13 A. Moreover, the possible values 
of [ug, r) are in the upper- right quadrant as shown in fig. 1241 

We see at the light of the three years WMAP data ref. |5j that hybrid inflation in the blue tilted regime rig > 1 is 
ruled out [ea. l4.2|l ] Q. That is, hybrid inflation in the regime jj? < /icrit i^ ruled out as well as hybrid inflation in the 
regime A^ > ylnt with x(0) > x(0)cr»t- 

Hybrid inflation in the red tilted regime /i^ > A^rit and x(0) < x{Q)crit fulflUs the three years WMAP value for ng 
ea. (l4.2|) . as well as the bound on the ratio r ea. H4.3|l . We can read from flg. l2Ulthat 

0.2 > r < 0.14 for 0.952 <ng< 0.97 . 

In addition, we flnd in flg. 1221 negative values for the running in this range, that is: 

-O.OOK dn^/dln/c < for 0.952 < n^ < 0.97 . 

There are clearly two regions which are not covered by hybrid inflation with only one inflaton field, neither by new 
inflation as shown in fig. |^ 
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FIG. 17: Hybrid inflation. The index n^ vs. A for ^^ = 1.7 A > ^^rit, 9^ = \ and x(0)crit ~ 5.7 \/A < x(0) < 10.7 /A. Notice 
that here we have both n^ > 1 and ris < 1 depending on the values of A and x(0)- ^U curves reach asymptoticaUy n^ = 1 for 

A -^ oo. 
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FIG. 18: Hybrid inflation. The ratio r vs. A for [j? = 1.7 A > ^^rit, 3^ = i and 2.7 \/A < x(0) < 5.7 VA ~ x(0)c 
decreases monotonically with A in this regime x(0) < x(0)crit and asymptotically vanishes for A -^ oo. 
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FIG. 19: Hybrid inflation. The ratio r vs. A for ^^ = 1.7 A > ^^^^t, g^ = j and x(0)crit ~ 5.7 \/A < x(0) < 10.7 ^/X. r 
shows an osciUatory pattern in this regime x(0) > x(0)crit and asymptotically vanishes for A ^ oo with no oscillations. The 
oscillations show up and concentrate with increasing amplitude for small A. 
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FIG. 20: Hybrid inflation. The ratio r vs. n^ for ^l^ = 1.7 A > ^^^^t, 5^ = ^ and 2.7 v^A < x(0) < 5.7 \/A ~ x(0)cr«t. Notice 
that here rig < 1 for all values of A and this range of x(0). We see that 0.2 > r > 0.14 for the interval 0.952 < n^ < 0.97. 
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FIG. 21; Hybrid inflation. The ratio r vs. Us for ^^ = 1.7 A > ^^^^t, g^ = j and x(0)crit ~ 5.7 \/A < x(0) < 13.7 \/A. Notice 
that we have here both n^ > 1 and n, < 1 depending on the values of A and x{0). AU curves end [x{0) < xiO)^,.^] or start 
[x(0) > xi^)crit\ at "-s = 1 which is the fixed point for aU values of x(0)- 
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FIG. 22: Hybrid inflation. The running dus/dlnk vs. n^ for /i^ = 1.7 A > ^^rit, S^ = i and 1.9 \/A < x(0) < 5.7 \/A ~ 
x(0)crit- Notice that we have here both positive and negative running dris/dlnk depending on the values of A and x(0)- Most 
values of the running are negative in this regime. For 0.952 < Us < 0.97, we have —0.001 < dUs/dlnk < 0. 
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FIG. 23: Hybrid inflation. The running dns/d\ak vs. n^ for ^^ = 1.7 A > fiirit, 9^ ^ \ and x(0)crit ^ 5.7 \/A < x(0) < 
13.7 VA. Notice that we have here both positive and negative running dns/d\nk depending on the values of A and x(0)- 
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FIG. 24: Regions described in the (r, n3)-plane by new and single field hybrid inflation for n^ < 1. The hybrid inflation border 
corresponds to /i^ — 1.7 A. For n^ > 1, all values of (r, n^) can be described by hybrid infiation (at least for r < 0.2, n^ < 1.15). 
The excluded region cannot be described by single field infiation (neither hybrid infiation, nor new infiation). Two or more 
fields infiation could describe such regions. 
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Simple single field inflation models have been recently studied within a numerical approach M^. In ref. [l^ 
cosmological data are fitted with the help of a Markov Chain Monte Carlo analysis. 

There is an interplay between the bounds of neutrino masses and the sign of (n^ — 1). A non-zero neutrino mass 
decreases the power in the small scales (large wavenumbers k). The same happens if n^ becomes smaller than unit. 
Therefore, if Us < 1, the power spectrum permits more stringent tests of the neutrino masses jl4| . The effect of 
neutrino masses for small scales Us > 1 can be cancelled by a spectral index n^ > 1. 
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